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Liebe DUG-Freunde, 

Da ist er nun, der erste Newsletter, der nicht in 
gedruckter Form erscheint. Ich wiinsche ihm eine 
freundliche Aufnahme. Die Reaktionen unserer 
Mitglieder auf die Anderung war sehr unterschied- 
lich (Seite 28). Ich danke fur das Verstandnis, das 
Sie der Umstellung entgegengebracht haben. Ich 
hoffe, dass der Inhalt dieses DNL Sie daon iiber- 
zeuget, dass wir unseren Standard halten. Die freie 
Mitgliedschaft bietet auch die Chance, unseren 
Kreis zu erweitern. Bitte weisen Sie CAS- 
Interessenten auf unsere Existenz hin. Alle Unter- 
lagen konnen auf der homepage - geanderte Adres- 
se im Editorial - nachgelesen werden. 

Falls es Probleme mit den pdf-Dateien geben sollte 
(Zeichensatze, ...), dann schreiben Sie bitte. Gleich- 
zeitig mit dem DNL#53 werden aktualisierte Aus- 
gaben der DNLs #1 und #2 aus 1991 zum Down- 
load angeboten. Auch die zugehorigen mth-Dateien 
wurden iiberarbeitet. Die Endungen v5, bzw. v6 
weisen auf Derive Versionen 5 und 6 hin. 

Es war reizvoll, aber auch sehr muhsam, die alten 
Dateien mit der modernen Textverarbeitung und 
Derive- Version wieder zum Leben zu erwecken. 
Ich habe mich bemiiht, die originale Erscheinungs- 
form so weit wie moglich unverandert zu belassen 
und dabei gleichzeitig die Inhalte zu aktualisieren. 
Machen Sie einen Blick hinein und freuen Sie sich 
iiber den Fortschritt des CAS. 

Ich wurde mehrfach ermahnt, dass ich in meinem 
Riickblick auf die ersten 13 Jahre DUG vergessen 
habe, mich auch bei meiner Frau Noor fur ihren 
Einsatz und ihre Einsicht fur die vielen von mir 
hinter dem PC verbrachten Stunden zu bedanken. 
Das sei hiermit "offiziell" nachgeholt. Ich mochte 
die Gelegenheit auch niitzen, meinen Dank an 
Bernhard Kutzler nachzutragen, er erwies sich 
immer als ein sehr hilfreicher Freund und Ratgeber. 

Beachten Sie bitte das reichhaltige Programm der 
Derive & TI-CAS, bzw. ACDCA-Konferenz 
TIME-2004 in Montreal. Im Rahmen dieser Konfe- 
renz wird ein Derive User Group Treffen stattfm- 
den. Es ware schon, dort viele Freunde zu treffen. 


Dear DUG-Community, 

Here it is, the first newsletter which does not ap- 
pear in printed form. I hope that it will be wel- 
comed by you. The reactions of our members on 
the change of publishing were very different (see 
page 28). I am very grateful for your understand- 
ing. The free membership gives the chance to find 
new members and contributors. Please inform 
CAS-interested people about the existence of the 
DERIVE & CAS-TI -Usergroup. All information 
can be downloaded from the homepages (one ad- 
dress has changed - see in Editorial). 

If there are problems with the pdf-files (type fonts, 
....) then please contact me. Together with DNL#53 
I offer adapted issues of DNL# 1 and #2 from 1991 
for downloading. The respective MTH-files were 
also revised and can now be used with Derive 5 and 
Derive 6. I hope that the contents of DNL#53 will 
convince you that we tried to keep our standard. 

It was a rewarding task, but also a tough one to 
revive the old files using modem text processors 
and recent Derive versions. I tried to leave the 
original appearance as unchanged as possible and 
simultaneously to revise and updroupate the con- 
tents for compatibility with Derive of today. Have a 
look and enjoy the evolution. 

Many of you reminded me that in my view back to 
the past I forgot to thank my wife Noor for her 
work and for her understanding for my many hours 
spent in front of the PC. Fd like to make this good. 
And Fd like to take the occasion to express my 
gratefulness to Bernhard Kutzler, he always was a 
very helpful friend and adviser. 

Please notice the rich program of TIME-2004 in 
Montreal (Derive & CAS-TI and ACDCA- 
Conference). We will have a DERIVE Usergroup 
meeting in the frame of this conference. It would be 
great to meet there many of you. 

Finally Fd like to call your attention to our friend 
David Halprin's idea to form a group of interested 
people to spread the work of the mathematician 
Ernesto Cesaro (page 44). 


Ich mochte Sie auch auf den Wunsch unseres 
Freundes David Halprin hinweisen, der mit 
einer Gmppe von Interessierten dem Lebens- 
werk von Ernesto Cesaro die gebuhrende 
Aufmerksamkeit widmen mochte (Seite 44). 

Mit herzlichen GriiBen 


With my best regards as ever 



DUG-Meeting 2004 
TIME-2004, Montreal 

Saturday, 17 July 2004, Lunchtime (tentative schedule) 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the TI-89/92/Voyage 200 as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS -TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: June 2004 

Deadline 15 May 2004 


Preview: Contributions waiting to be published 

Finite continued fractions St. Welke, GER 

Some simulations of Random Experiments, J. Bohm, AUT 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encri ption, J. Bohm 

CAD-Design with DERIVE and the Tl, J. Bohm 

Avoiding Convolution and Transforming Methods, M. Lesmes-Acosta, COL 
Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Fuzzy Logic, G. Hagen, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 
Modelling real date: enthalpie values 
A TVM-Calculator for DERIVE, M. R. Phillips, USA 
Shares, Put and Call, J. Bohm, AUT 

and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT, 

Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung : F achzeitschrift 
Herausgeber: Mag Josef Bohm 
Herstellung: Selbstverlag 

URLs: http : / / www . austromath . at/ dug/ index . htm 

http : //shop . bk- teachware . com/k . asp?session=2 57 9115&kat=ll 
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Manuel Garcia Mateos 


garcia_mateos@WEB . DE 


Limit of a Sequence 

Dear Derivians, 

I need to know how to calculate a limit on a recursive sequence! For example I defined the following 
sequence mittelwert (n) : =iterates ( [b, (a+b) /2 ] , [a, b] , [0,1], n) 

I received a list and I know that the sequence mittelwert (n) converges to 2/3. 

I think there is no problem to calculate the limit of a explicit sequence, but what about a recursive? 


Thanks a lot for your help! 
Manuel 


Stefan Welke Spwelke@AOL.COM 

Dear Manuel, 

you can obtain your sequence of vectors by matrix-multiplication: 

f_0_ ( [a, b] ) = [[ 0 , 1 ],[ 1/2 , 1/2 ] ] * [a, b] and generally: 
f _n_ ( [ a , b ] ) = [ [ 0 , 1 ] , [ 1/2 , 1/2 ]] A n*[a,b] 

You can compute the eigenvalues l_1 , l_2 and eigenvectors e_1 , e_2 of this matrix and similar matri- 
ces, which describe recursive sequences. Then you can deompose [a,b] as a linear combination 
[a, b] = alpha*e_l + beta*e_2 and by linear algebra you obtain: 

f _n ( [ a , b ] ) = alpha * l_l A n * e_l + beta * l_2 A n * e_2 , 

which is the sum of two geometric sequences. 

I hope this helps. 

Stefan Welke 


Ignacio Larrosa Canestro ilarrosa@MUNDO-R.COM 

a(n+l) = (a(n)+ a(n-l)/2 ===> 2a(n+l) = a(n) + a(n-l) 

2a (n+1 ) + a(n) = 2a (n) + a (n-1 ) , for each n >= 1 
Then, 

2a(n+l) + a(n) = 2a(2) + a(l) 

As a(n+1) is between a(n) and a(n-1), and |a(n+1) - a(n)| = |(a(n) - a(n-1)|/2, a(n) has limit as n — > inf. 
Let it a. Then 

2a + a = 3a = 2a(2) + a(l) ====> a = (2a(2) + a(l))/3 

Saludos, 

Ignacio Larrosa Canestro 
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Wim De Jong wimdeJongTwo@A0L.C0M 

Hi Manuel, 

You are dealing with the sequence of vectors (a(n),b(n)) defined by 


(a(n+1 ),b(n+1 ))=(b(n),0.5(a(n)+b(n))) with (a(0),b(0))=(01). 

So you have a(n+1)=b(n) and b(n+1)=0.5(a(n)+b(n)), i.e. 2b(n+1)-b(n)-b(n-1)=0. This homogeneous 
linear recurrence relation (with constant coefficients) of degree 2 is treated in the same way as a ho- 
mogeneous 2nd order linear differential equation (w.c.c.). Try b(n)=t A n. This yields the auxiliary equa- 
tion 2t A 2-t-1=0, which has the solutions -0.5 and 1. So the general solution of the recurrence relation 
is b(n)=A(-0.5) A n+B, where A and B are arbitrary constants. Hence (a(n),b(n)) = (A(-0.5) A (n-1)+B, 
A(-0.5) A n+B). Now (a(0),b(0))=(0,1) yields the special solution ((1/3)(-0.5) A (n-1)+2/3,(1/3)(- 
0.5) A n+2/3). And as n tends to infinity (a(n),b(n)) tends to (2/3, 2/3). 


Another interesting technique for solving 2b(n+1)-b(n)-b(n-1)=0 with b(0)=1 and b(1)=0.5 is the gen- 
erating function method. Consider the so-called generating function G(x)=b(0)+b(1)x+b(2)x A 2+ 
b(3)x A 3+ ... (this is not really a function; it is an expression for b(n) which allows convenient algebraic 
operations). 


G(x) = 1+0.5x+0.5(b(1)+b(0))x A 2+0.5(b(2)+b(1))x A 3+ ...= 

= 1 +0.5x{1 +b(1 )x+b(2)x A 2+...} + 0.5x A 2{b(0)+b(1)x+b(2)x A 2+...} = 1+0.5xG(x)+0.5x A 2G(x), 


which yields 

G(x)=1/(1-0.5x-0.5x A 2)=1/{(1+0.5x)(1-x)} = (use partial fractions) (1/3)/(1+0.5x)+(2/3)/(1-x) = (use the 
binomial expansion) (1/3){1-0.5x+(0.5 A 2)x A 2-(0.5 A 3)x A 3+...}+(2/3){1+x+x A 2+...}. In this formal power 
series the coefficient of x A n (n=0,1, 2, ...) is (1/3)(-0.5) A n +2/3. 


Thus b(n)=(1/3)(-0.5) A n+2/3 for n=0,1 ,2 

Cheers, Wim 


Patrick West patrick@fpwest.com, www.fpwest.com 

Folks, 

I just recently moved to Pocket PC from the HP 200LX palmtop. 

I am happy to report that my old DOS version of Derive works nicely under PocketDOS on my 
Ipaq 3845. 

I also obtained the ti-89 emulator shell for Pocket PC and loaded the rom from my 89 and it works 
nicely. 

James F. Gordon jfgordon@FRONTIERNET.NET 

Is Derive 6 able to perform implicit differentiation? For example, how can I find dy/dx for the following 
implicit function of x: 

y A 3*sin(x*y) = x A 4.5 

I was not able to find implicit differentiation in the introductory text that accompanied the CD. 

Thank you for any help you can give me on this. 

James 

Louis F. Lowell 

Does anyone know how to do implicit differentiation with Derive? Specially something like 
x + sqrt(x*y) - y = 4 or y A 2 + e A (x*y) = sin(x A 2*y A 2)+2 
Thanks for any help anyone can provide. 
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Aleksey D. Tetvorko chib@MEGASTYLE.COM 

I cannot do best then you can see in implicit .mth. 

Aleksey 


===============implicit .mth : 

y A 3*SIN (x*y)=x A (4.5) 


DIF (y A 3*SIN (x*y) =x A (4.5) , x) *t+DIF (y A 3*SIN (x*y) =x A (4.5) , y) *u 
; Simp (#2) 

(u*x*y A 3+t*y A 4) *C0S (x*y) +3*u*y A 2*SIN (x*y)=9*t*x A (7/2) /2 
; Solve ( #3 , u) 

SOLVE ( (u*x*y A 3+t*y A 4) *C0S (x*y) +3*u*y A 2*SIN (x*y)=9*t*x A (7/2) /2 , u) 
; Simp (Solve ( #3 , u) ) 

u=t* ( 9*x A (7/2 ) -2*y A 4*COS (x*y) ) / (2*y A 2* (x*y*C0S (x*y) +3*SIN (x*y) ) ) 
u/t 


SUBST (u/t,u,RHS (u=t* (9*x A (7/2) - 

2*y A 4*c0S (x*y) ) / (2*y A 2* (x*y*C0S (x*y) +3*SIN (x*y) ) ) ) ) 

; Simp (#7) 

( 9*x A (7/2) -2*y A 4*COS (x*y) ) / (2*y A 2* (x*y*C0S (x*y) +3*SIN (x*y) ) ) 
===============implicit . mth================================= 


MacDonald Phillips phillipsm@GAO . gov 

There is a command in Derive 6 for implicit differentiation. It is in one of the utility files and will be 
automatically called up when you type it in. First, bring the x A 4.5 to the left hand side so you have an 
expression that is implicitly equal to zero. The command is: 

IMP_DIF(u,x,y,n) 

This gives the nth implicit derivative of y wrt to x of expression u. Since n defaults to 1 all you have to 
enter for your example is: 


Hope this helps. 


IMP_DIF (y A 3*sin (x*y) -x A 4 . 5, x, y) 


Terence Etchells 

imp dif (u, x, y, n) can be used also in earlier Derive versions (function of DIFF_APPS.mth 
utility file) 

IMP_DIF(x + 4(xy) - y = 4, x, y) = ? 


2 xy 2 2 

IMP_DIF(y + e = SIN(x + y ) + 2, x, y) = 


xy 

ye x 


xy y 

xe + 2y 


It's obvious , that it doesn't make any sense to apply imp_difon equations. Students must know this 
and bring all expressions on one side of the equation (according to Don's mail) and are not allowed to 
believe the result without any doubts. In case if there is a constant on one side of the equation Derive 
5 returns a ? showing that something must be wrong. Derive 6 does not and offers a "result". 


IMP_DIFCx + VC x-y) 


y = 4) = 


' y 

^ X 




-I 
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It is very interesting to follow Stepwise Simplification of Implicit Differentiation of DERIVE 6. 
See the first step: 

3 4.5 

#1: IMP_DIF(y -SIN(x-y) - x ) 

d d d 

— (F(x) + G(x)) * — F(x) + — G(x) 

dx dx dx 



fd ^ 

fd 9/2 d 3 

d 

9/2 d 3 



— dk - 

— C- x ) + — (y -SIN(x-y)) 

i 

/-N 

1 

X 

) + — (y -SIN(x.y)) 


d 

^dy > 

^dx dx > 

dx 

dx 



#2 : ITERATE 


dk 


dx 


d 3 9/2 

— (y -SIN(x-y) - x ) 
dy 


dk. 


d 3 9/2 

— (y -SIN(x.y) - x ) 
dy 


and the final ones: 


#17: ITERATE 


7/2 

9*x 4 

y *C03(x*y) 


y *COS(x*y)*x + 


( d 3 

— y 

.dy 


0 


*SIN(x*y) 


d n n - 1 

— x =& n - x 

dx 


#18: 


#19: 


ITERATE 


0, 


d*k , 


r 


ITERATE 


0, 


d*k , 


7/2 \ 

9*x 4 

y *COS(x*y) 

2 

, 0 

3 2 

y *COS(x*y)*x + 3 *y *SIN(x*y) j 

7/2 4 

9*x - 2-y *COS(x*y) 

, 0 

2 

2*y * (x*y*COS(x*y) + 3*SIN(x*y)) 


I tried to extend imp_dif() to imp_diff(), which works also for implicitly given expressions, without col- 
lecting both sides on on side of the equation: Josef 


imp_diff(u, x, y, n := i) := 

If RHS(u) = 0 

ttl : IMP_DIF(u, x j. y, n) 

IMP_DIF(LHS(u) - RHS(u), x„ y„ n) 
IMP_DIF(LHS(u) - RHS(u) „ x„ y„ n) 

y (sT(x y) + 2 x) 

#2: imp_diff(x + T(xy} - y = 4) = 

x (2 y - sT(x y>> 

y (T(x y) + 2 x) 

#3: inp_diff(4 - x - T(xy) + y) = 

x (2 y - T{x y)) 


The following is a letter from Kathrin who needed support in programming: 

Hello all! 

For my skilled work I have to program a function (in Derive) which computes the error of Taylor series. 
This function which I called "f ehlerfunktion" works: 
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fehlerfunkt ion(funktion , n, epunkt, F) := 

Prog 

F := funktion - TAVLOR(funktion , x, epunkt, n) 
RETURN F 


(funktion is e.g. sin(x) and epunkt is the the position at which the polynomial shall approximate) 


Then I tried to compute the error at the point bpunkt because I am not allowed to set x : =bpunkt 
and to use the function solve in a second step. I had to write a new function which I called fehler 
and this one does not work and I don't know why: 

FEHLER(funktion , n, epunkt, bpunkt, M) := 

Prog 

M := funktion - TAVLOR(funktion , x, epunkt, n) 
x := bpunkt 
RETURN M 

My third job was to find a function which computes the n which approximates at bpunkt as good as 
the user wants (the acceptable error I called fehlerwert): 

napproximat ion(funktion , n, epunkt, bpunkt, fehlerwert) := 

Prog 

Loop 

If FEHLER(funktion , n, epunkt, bpunkt) > fehlerwert 
RETURN n 
n :+ 1 

Because my second function does not work I don't know if the third one is correct (using this function 
Derive just computes and computes without displaying a result) . 

This is my first time using Derive and maybe my questions sound very banal. Nevertheless I would be 
very happy if you could help me! Thanks for any kind of help anybody can provide! 

with kind regards 
Kathrin Becker 


DNL: 

Dear Kathrin, 

I wanted to follow your ideas and tried to correct your programming efforts: 

FEHLER ( funktion, n, epunkt , bpunkt) : = 

SUBST ( funktion -TAYLOR ( funktion, x, epunkt, n) , x, bpunkt) 

and then using this auxiliary function use the LOOP as follows: 

napproximation(funktion, epunkt, bpunkt, fehlerwert, n := G) := 

Loop 

If ABS(FEHLER(funktion , n, epunkt, bpunkt)) < fehlerwert 
RETURN n 
n :+ 1 

which must be typed in as: 

napproximat ion (funktion, epunkt, bpunkt, fehlerwert, n : =0 ) : = 

LOOP (IF (ABS (FEHLER (funktion, n, epunkt, bpunkt)) < fehlerwert, RETURN n) , 
n :+ 1) 


Then it works as you can see below. Johann Wiesenbauer provided a very compact function, which is 
more general because it allows also variables t x. 

Josef 







Johann Wiesenbauer 

Hi Kathrin, 

For a start, try the following routine 

napprox(u, x, xO, xl, epsilon, n := 0) := 

LOOP( 

IF(ABS(SUBST(u - TAYLOR(u, x, xO, n), x, xl)) < epsilon, RETURN n), 
n :+ 1) 

and tell me, if it does what you wanted it to do. 


Cheers, Johann 


1 

napprox SIN(x) , x, 1, 2, = 6 

l 1000 . 


napproximationfSIN(x) , 1, 2, 0.0001) = 6 


Some weeks ago I received an e-mail from Kim Hendrickx, Tl-Support Centre, who forwarded a prob- 
lem, which was sent to him recently: This was the problem: 

Prof. Dr. P. 

Prof P. wants to plot the graph of 

y = ln(cos(x)) / ln(cos(x)) over the domain Db(f) = {x | -pi/2 + 2 k* pi <x< pi/2 + 2k* pi , k integer). 

The function value is 1 where the function is defined. 


On the TI-83Plus: 



The graph is correct! But .... 


On the TI-89, TI-92Plus, Voyage 200: 



The CAS-Calculator plots the line y = 1 and represents incorrectly function values outside of the do- 
main. But if I change the function definition slightly then: 
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This was my answer and I would like to ask all of you about your opinion , because I am not sure if I 
am right with my conclusions. Josef 

Dear Professor P., Dear Kim, 

Many thanks for the interesting graphics problem, which seems to be more a mathematical than a 
graphical problem. 

First of all, DERIVE and MuPad are both plotting;; = 1 as graph of your function. 

I believe that there are two reasons which are responsible for the graphic representation of this quo- 
tient-function of two logs. (By the way, ln(x)/ln(x) shows also y = 1, which should have a graph only 
for x > 0 following your intentions??) 

(1) The Computeralgebra System seems to recognize numerator and denominator as identical and can- 
cels before evaluating numerator and denominator .... 

(2) More likely the CASs are calculating with complex numbers. 

Before demonstrating a possible solution of your problem, i.e. to reproduce the TI-83 plot on the 
TI-CAS-Calculators Ed like to mention some remarkable results, which could inspire discussion in 
classroom. 

Investigating the "wrong” graph using TRACE you will find that there is no function value given for 

all places x with cos(x) = 1 (log = 0). This happens for x = 0, 2n , 4n, but for x = n, 3tt, .... we find 

non-real result. If we change Complex Format (MODE) from REAL to RECTANGULAR, we 
obtain again the quotient 1, because of ln(-l) = n*i. 

The TI-83 cannot perform all these calculations, because of its lack of computeralgebra. I believe that 
this calculator tries to evaluate numerator and denominator numerically and being unable to do this in 
the complex region leaves the function value undefined - and consequently doesn't plot any pixel. 


I tried to produce the desired plot and wrote a function which prevents cancellation and and simulates 
the TI-83+ more or less. Function Incos(x) is written into the Y=-Editor and then plotted. 



Complex Format REAL Complex Format: RECTANGULAR 


Playing around with this function I came across another interesting fact. Increasing the xres-value in 
WINDOW, the resulting plot shows a gap (see the picture below). But this happens only for xres = 5 
with the ZoomTrig-setting. Exchanging cos by sin, the gap appears in the left halfplane. 

Any explanation?? 

^qqJt raceyRegraphyMaihlpraLolT yWt I hope that my attempt to interpret this interesting phenomenon 

could help and remain .... 


I MAIN 


FiliD flUTD 


rune 


Unfortunately I only received an very friendly answer from Kim, 
but non from Prof P. So I don't know if I could satisfy him. 
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Gerhard Pachler, Austria 

The students solved the problem to rewrite the following fraction without roots in the denominator: 

30 


Calculation by hands reults in: 


3 V 5 + 5 V 3 + 2 V 30 

Then we repeat the calculation using the TI-92: 


3 V 5 + 5 V 3 + 2V30 

30 =4i + Jl-2sT2 


f^JShuTbnalc^IotrenlPn^Iofcie^n UpT 1 

■ 3 JB + 5 ■ J3 + a 

3-J5 + 5-J3 

■ 2-J30 ■* b 

2 J30 

. 30 

30 

a + b 

[2-j6 + 3)-j5 + 5-j3 

30 (a - b) 
Ca + b) (a - b) 

J5+J3-2J2 

30*<a-b>/<Ca+b>*ta-b;>>l 

MAIN RAD HLITD 

FUNC 1/30 


I'FiT^riT FF^ Y F2 t Y FHt Y F? T Ffi^ T "l 
Rlgebra|Calc|Oiher|PrgNlO|Clean Up| J 


50 (a - b) 

(a + b) (a - b) 

50 _ 50 (a ~ b) 


12-J6 + 3J-J5 + 5-J3 
J 5 +J 3 - 2 J 2 


a + b (a + b) (a - b) 
■ + J3 - 2 -J2 


30 


(2 -JO + 5)-j5 + 5--J5 


false 

false 


...=30/<X2*J<6 >+3>*J <5>+5*J<3>>l 


Why does the Tl not recognize the equivalence of the expressions? How can I convince the machine? 


DNU 

This is my treatment of the roots (fortunately not at the dentist!) 

After some tries which were answered with false / fac- 
tored out sqrt(15) in numerator and denominator and 
cancelled the root - and now the equivalence became 

true// 

Then I proceed: 

Copying the two expressions as ans(2) and ans(1) into 
the edit line and then approximating ans(2) - ans(1) / 
receive a result which differs from 0 - which might be an 
internal numerical affair?? 

But if I copy the expressions into the edit line or if I type 
them in the Exact Mode refuses to perform any calcula- 
tion, approximating the difference I obtain 0.00000.... for 
all decimal places. 

Evaluating the equation exactly I receive false, but ap- 
proximating it I get true again. 

So I am not surprised that approximating ans(2)=ans(1) 

is false. 

It seems that we could write a PhD on the behaviour ot 
the CAS-Tls in treating expressions containing roots. 


Best regards 
Josef 



: : Kfcar;.|prS,I0lc:,:S *1 1 

' 

■ 

■ 

■ 

2-JTB 

J5 +f3 + 2J2 
30 

3-J5 + 5-J3 + 
J5 + J3 - 2-J2 

J5 +f3 - 2 f2 

= j5+f3-2j2 true 

r— 1 . 13969166032 

2 ■ -J 30 

1 . 13969166032 

ans<2>-ans<l> 

MAIN RAD EXACT FUNC 10/10 



^•^rgniofiuiir Uf-T | 

■J5 +f3 - 2-f2 

1 . 13969166032 

. 30 

30 

3-J5 + 5-J3 + 2 J30 

■ J5 + J3 - 2-J2 
. 30 

(2-J1 + 3) j5 + 5 J3 
JB + J3 - 2 j2 
-(J5 + J3 - 2-J2) 

(2f6 + 3)f5 + 5j3 



ans<2>-ans<l> 

MAIN RAD EXACT 

FUNC H/10 


f?^hL^r,kVkh^r1Pr^I0l 


□ 

— i . uijud 

■ _ (f5+j3 2 

3-J5 + 5-J3 + 2-J30 K 

0. 

■ j = ^ j= = JB + J3 - 2 

3-J5 + 5-J3 + 2-J30 

■ f= ^ j= = JB + J3 - 2 

3f5 + 5f3 + 2-J30 

innnnnnn t TTj 

.■m 

00000000000 
■J2 false 

J2 EHfE 

ans<2>-ans<l> 

MAIN RAD EXACT FUNC 1/10 


Are there any other findings?? 
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Mate Matica 

The expression (see page 3) abs(int(#e A (#i*pi*x A 2)*#e A (#i*2*pi*f*x),x,-10,10)) equals 


ERF 

' 42-4n(1 - i) ■ (f - 10) 

- ERF 

' 42-4n-(1 - i) ■ (f + 10) ] 

2 

2 




2 


(it can be obtained using Derive with some tricks). Now the numerical errors (generated by I NT) are 
eliminated and the plot is ok. 



Volker Loose 

Hello all, 

is there a function in Derive which returns the length of the periodicalpart of a decimal fraction, 
is there a function in Derive which returns the length of the preperiodical part of a decimal fraction? 

Example: 1/28=0.03571428571428... 

preperiodical part 03, length 2 
periodical part 571428, length 6 

Thanks 

Volker 


Johann Wiesenbauer 

Actually this is the same question Peter Schofield asked me some time ago. 

I promised to write a utility function then, but sadly enough forgot about it. Sorry, Peter! Well, at last in 
the attached file there is a first version which might live up to the expectations of both of you. 

Cheers, 

Johann 


Albert Rich 

Hello Johann, 

In your Titbits 26 in DNL#52 you challenged the reader to determine the number of digits in the re- 
cently discoved largest known Mersenne prime: 2 A 20996011-1 using Derive. That's easy, simply 
approximate it giving 

1. 2597 68954- 10 A (6. 32042 9-10 A 6) 

The number of digits is this power of 1 0 plus 1 . 

Probably the easiest challenge you have ever issued. :=) 

Aloha, Albert 
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On Josephus’ Permutations 

Stefan Welke, Bonn, Germany, spwelke@aol.com 


1 Faster computation 

The following code presents two functions, which compute Josephus’ permutaions. These permuta- 
tions were described in DNL #52 by Riideger Baumann as a programming challenge. Riideger’ s func- 
tion is named josephusl(n , k). You can find the code in DNL #52. 


josephus2(n, k) = 

Prog 

[v := [1, ■ ■ ■ , n], w := [], q := n, p := 1, j := 1] 

Loop 

#1: If q = 0 

RETURN w 
If p = k 

[w := APPEND (w, [v|j]), p := I, v := DELETE (v , j), IF(j = q , j := I) p q 1] 
[p :+ 1, IF(j < q, j :+ i, j := 1)] 


j osephus2 (n, k) is slightly faster than Riideger’ s function, but it shows the same time behaviour 
with respect to the arguments n and k : The computation time depends nearly linearly on n and on 
k . 

josephus3(n, k, v, w, q, p, j, m, z) := 

Prog 

[v := [1 n], w := [], q := n, p := 1, j := 0, in, z, v] 

Loop 

If k > q exit 
m := FLOOR ((q - j)/k) 
z := VECTOR (j + i*k, i, I, m) 
w = APPEND(w, viz) 
v := DELETE (v , z) 
j := m*k + j - q 

q := q - m 

j := I + MODCj + k - I, q) 

Loop 

If q = 0 
RETURN w 

w := APPENDCw, [vtj]) 
v := DELETE (v , j) 

q - 1 

j := I + MODCj - 2 + k, q) 

The second version is reasonably faster than the preceeding one. The time behaviour depends mainly 
on n and not on k . The following table shows some computation times in seconds on my 3.6 Ghz 
Pentium. 


n=300 k=20 k=40 k=150 k=200 k=280 k=400 k=600 


Derive 6 provides a new DELETE-command: 
DELETE ( [a,b, c, d, e] , [1,3,5]) removes 
elements a, c and e from the first vector. 


josephusl 1.77 3.55 13.3 17.4 24.5 34.9 52.4 
j osephus2 0.250 0.453 1.58 2.09 2.91 4.14 6.28 
j osephus3 0.0150.0310.0470.0470.0470.0470. 047 
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2 The set of special Josepus permutations 

The Josephus functions generate lists of numbers, which represent permutations of the set {1,2,..., 72 } 
of the first n natural numbers. Now the question is, if every permutation is a Josephus permutation. 
The number of different pemutations of n different objects is n\ , so the sequence 
{josephus(n,l), josephus(n,2),..., josephus(n,\ + n\)} must contain at least two equal permu- 
taions, and the next example demonstrates this fact: 

#7: josephus3(20 , 1+20!) 

#3: [1, 2, 3, 4, 5, 6, 7, 3, 9, 10, 11, 12, 13, 14, 15, 16, 17, 13, 19, 20] 

But the following theorem shows, that Josephus permutations are only a subset of the group of all 
permutations for n > 4 , because lcm(l,2,..., n) <n\ for n > 4 . 

Theorem 

The equation josephus{n , k ) = josephus(n , k + lcm(l, 2, n )) holds for all n,k e N* . 

#9: lcm ( [ i , . .., 200]) 

#10: 3372935888329262 64639465766794841407432394382785157234228847021917234018060677390066992000 


#11: josephus3(200 , 1 + LCM([1, 200])) 


#12: 

[1, 2, 

3, 4, 

5, 6 

, 7, 

8, 9, 

10 , 

11, 12 

• 1 13 

, 14, 

15, 

16, 17 

, 18, 

19, 

20, 

21, 22 

, 23, 

24, 

25, 26, 27, 

, 28, 


163, 

164, 

165, 

166, 

167, 

168, 

169, 

170, 

171, 

172, 

173, 

174, 

175, 

176, 

177, 

178, 

179, 

180, 181, 

182, 


183, 

184, 

185, 

186, 

187, 

188, 

189, 

190, 

191, 

192, 

193, 

194, 

195, 

196, 

197, 

198, 

199, 

200] 



(Computation time: 0.031 sec) 


3 Encryption 

A slight modification of j osephus3 (n, k) allows us to permute the order of characters in a given 
string by a Josephus permutation, where n is the number of characters in that string. 

josephus4(1ist , k , v, w, j, q, in, z) = 

Prog 

[v = list, w := [], j = 0, q = DIM(list)] 

Loop 

If k > q exit 

in = FLOOR (Cq - j)/k) 

z = VECTOR ( j + i-k, i, 1, in) 

w = APPEND (w, viz) 

v = DELETE (v , z) 

#13 : j := m*k + j - q 

q = q - m 

j = 1 + MODCj + k - 1, q) 

Loop 

If q = 0 
RETURN w 

w = APPENDCw, [vij]) 
v = DELETE (v , j) 

q - J 

j = 1 + MOD( j - 2 + k, q) 

#14: josephus4([a , b, c, d, e], 5) 


#15: 


[e, a, c, d, b] 
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We now use this function to define encoding and decoding functions for strings: 


#16: encodeCstring, k) = C0DES_T0_NAME(josephus4(NAME_T0_C0DES( string) , k)) 

#17: encodeCHallo, fellow derivians, 11) 

#13: lnlswlrfavdoe iH o,aeil 

decodefstring , k, w, v, txt) := 

Prog 

#19: w = NAME_TO_CODES( string) 

v = josephus3(DIM(w) , k) 

COD E S_TO_N AM E ( SORT ( [ v , w] 1 ) 1 i2) 

#20: decode(lnlswlrfavdoe iHo,aeil, 11) 

#21: Hallo, fellow derivians 

Permutations are no good means for encryption. But here the whole text is permuted, and not only the 
alphabet. Thus the longer the text, the more possible permutations are at hand. For example: A text of 
200 characters admits LCM([1,...,200]), which is about 3. 372935888- 1 0 89 , possible Josephus permuta- 
tions. 

The given decoding function works painfully slow, due to the two matrix tranpositions in the code. 
Function decode2 (string, k) works considerably faster: 


decode2 (string, k, v_, n_, q_, z_, w_) = 

Prog 

v_ = NAME_TO_CODES(string) 
n_ = DIM(string) 

q_ := I 

Z_ := [I n_] 

#22: w_ = josephus3(n_, k) 

Loop 

If q_ = n_ + 1 

RETURN COD E S_TO_N AM E ( z_) 
z_ = REPLACE (v_iq_, z_, w_ iq_) 
q_ :+ 1 

#23: deoode2(lnlswlrfavdoe iHo,aeil, 11) 

#24: Hallo, fellow derivians 

Decoding a string is easy, if we know the number of the permutation given by the Josephus function. 
Here is my challenge for the readers of the Derive News Letter : 


Decode the following string, which was encoded with josephus4 and the resulting encoding 
function. 


"ruimunlasee, dvretol -uo noqudne avmdaaW eowmckadwao h,ot ni e ei aamo 
ntphnm sobayoc, o,p,or eeoi nn 1 stpe pogut rf tauot edstyiplwrdnl , n gpi 
dhnselt or as iib cnrnamln a .onlae f naha a rydi ra,ierg gent .naemnidgo 
yppg rmp m eongdrhoy yrpetOa pshun enseogadthorti ceo te h miryeryynrvl 
ii dtc leT lan f ddge iade-r.ermm p " 


There is no space character at the end of the first four lines. 
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Der Artikel von Franz SCHLOGLHOFER im DNL52 hat mir ausgezeichnet gefallen, 
ebenso Ihre Erganzungen. Es ist mir gelungen, die ’’threads” fur vier Punkte, also fur 
die Bezier-Kurve 3. Grades in der Ebene hinzubekommen. Wenn Sie Zeit und Lust ha- 
ben, werfen Sie einen Blick darauf. 

Allerdings beschleicht mich der Verdacht, dass das einfacher auch gehen konn- 
te/miisste. Das ’’riecht" formlich nach einer Rekursion, oder? 

Gerhard Hagen 


Bezier Curve of degree 3 - Threads for four points 

Gerhard Hagen, Althofen 

We start with an example: A(-5|0), B(0|0), C(8|0), D(9|9) and develop threads for groups of three 
points: 


[A := [-5, 9], B := [0, 
threads2(p1 _ , p2_ , p3_ , 


0] , C := [8 , 0] , D := [9, 9]] 


n) := UECT0R 


p1_ + ( P 2_ - P 1_), 

n 


i 

p2_ + ( p3_ 

n 


P2_) 


0 , n 


Two threads2 are composed using 4 points and resulting in help2. 

help2(p1_, p2_, p3_, p4_, n) := APPEND( threads2(p1 _ , p2_, p3_, n) , threads2(p2_ , p3_, p4_, n)) 
Simplfying help2 we note that some threads appear twice. help3 resolves this ’’problem”. 

help3(p1_, p2_ , p3_, p4_ , n) := DELETE(help2(p1 _ , p2_, p3_, p4_, n) , n + 1) 


help3( A , B , C, D , 4) 
help3( A , B , C, D , 20) 




This auxiliary construction does not osculate the curve except in the end points. 

3 2 2 3 

bez3(a_, b_, c_, d_, t) := (1 - t) a_ + 3(1 - t) tb_ + 3(1 - t)t c_ + t d, 

bez3( A , B, C, D , t) 



Final composition of threads3 from threads2 and help3 requires some ’’Try and Error Testing 
and Playing”: 

The groups must fit together and groups of three points must be chosen correctly. In the respective 
thread2 only each i-th segment is tangent of the curve. 
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threads3(p1 _ , p2_, p3_, pH_, n) := UECTOR 


( threads2( (help3(p1 _ , p2_, p3_, p4_, n)) , (help3(p1_, p2_ 

i , 1 


i , 2 


(help3(p1_, p2_ , p3_ , p4_, n)) 

i 

+ n,2 

n)) 

i 

, i , 

threads3( A , 

B, C, 

D, 

*0 

threads3( A , 

B, C, 

D, 

20) 



An open question: What is the locus of the endpoints of the threads? 

Albert Rich & Johann Wiesenbauer: Facing the Challenge 

Hello Johann, 

In your Titbits 26 in DNL#52 you challenged the reader to determine the number of digits in the re- 
cently discoved largest known Mersenne prime: 2 A 20996011-1 using Derive. That's easy, simply 
approximate it giving 

1 .259768954 • 10 A (6.320429* 10 A 6) 

The number of digits is this power of 1 0 plus 1 . 

Probably the easiest challenge you have ever issued. :=) 

Aloha, Albert 


Hello Albert, 

It's certainly very impressive that Derive can cope in approximate mode with the problem of determin- 
ing the number of digits of this huge number, in particular in the light of the recent discussion about 
APPROX on the newsgroup. I have a slight preference for my solution though, which I gave at the end 
of the my ff Titbis(26) ?f for arbitrary Mersenne numbers 2 A p-l, namely 

digits (p) :=floor(p*log(2,10))+l 

for two reasons: 

- This method will still work for even higher exponents without adjusting the precision (in other 
words, FLOOR will take care of that adjustment itself!) 

- It can be asily adapted to other bases b>2 (you only have to replace 10 by b in the formula above) 

Nevertheless, I agree that this challenge wasn't too hard for a change, in particular, if you can resort to 
such a powerful tool as Derive :)) 

Happy Easter, 

Johann 
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Jan Vermeylen, Belgium 
Hello Josef and Albert, 

Can I bother you with a problem I got from an engineer? 


He wants to estimate the value of the parameters K and k in an expression where he does experi- 
ments measuring x and t, while changing constants Ca, Cb, Cc, Cd for every experiment. 


The expression is: 


a = Ca(K- 1) 

b = ~(K ■ (Ca + Cb) + Cc+ Cd) 

Cc-Cd 


\i = K ■ Cb 
q 2 = b 2 


Ca 

4 -a-\i 


( 2 ax') 

... K , b-q 

with t = In — 

k q t | 2 ax 

l b + c l y 


I think this is a straightforward application of the FIT function, but to my surprise, I get no good an- 
swer. 

What do I do wrong? 


Jan successfully produced a expression for t depending on K and k using the Computeralgebra by 
substituting step by step. As you can imagine the resulting expression is a very bulky one. Finally Jan 
tried regression using given data from below. It didn't work - Derive Online Help: The de- 
pendence on the parametric variables should be linear. 



X 

ca 

cb 

cc 

cd 

expr 




1 

0.3 

0.5 

- 1.8 

2.3 

7.6 


Any idea, how to solve the problem?? 

data : = 

2 

- 1.9 

3.4 

2.4 

- 5.7 

- 0.21 




. 3 

-5.4 

7.4 

-0.3 

-0.45 

2.56 




A Macro for Derive 6 

Johann Wiesenbauer wrote a Macro for earlier Derive versions which happened to fail for upgraded 
Derive versions , because Derive became “more clever”. Johann used all his tricks to keep Derive not 
too clever for this special purpose. And again was a call from German teachers who wanted to have 
Johann's valueable tool adapted for using it with Derive 6. (See also DNL#31 and DNL#42.) 

Albert Rich's wrote 

In Derive 4 and 5, when an equation is simplified (as distinct from being solved), the left and right sides of the 
equation are simplified independently, and the results returned as an equation. So basically, the equality opera- 
tor is just a place holder. 

In Derive 6, when an equation is simplified, it is replaced by equivalent, but simpler, expression having the same 
Boolean truth-value. This is the way virtually all other expressions are simplified in Derive. In addition to 
being more consistent for the user, throughout the Derive 6 equation and inequality solver the assumption is 
made that equations will be fully simplified. 

Therefore, I wonder if it would be possible for you to modify your routine for step-by-step equation solving so 
that equations are represented as 2-element vectors. Stored in this way, the left and right sides of equations can 
be independently simplified and manipulated. 
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And it was possible for Johann: 


Johann Wiesenbauer 

Hello Derivers, 

On request of many users I have updated my macro for solving equations to DfW6 by taking into ac- 
count that version 6 handles equations a bit differently. In particular, I had to deal with each side of the 
equation separately. The output of every DO-command is now a 2-row matrix in order to suppress 
commas. I hope you can live with this. On the other hand, I have added a new command sqr(@) to 
compute square roots. In a similar way, other useful functions (e.g. for expanding and factoring poly- 
nomials) could be added, if there is a need for it. 

I wish you much fun with this macro and look forward to any comments or suggestions you might have 
as to further improvements! 

Cheers 

Johann 


Here is the file (it is very enjoyable to follow 
Johann's programming skills. Josef 


#1 : 


RECORD(u) 

Prog 

0 : = 

Pi : = 
P2 : = 

1 : = 
r : = 


□ 

LHS(u) 

RHS(u) 

[PI] 

[P2] 


u 


#3 : 


sqr(u, u ) := 

Prog 

If NUMBER? (u) 

#2 : RETURN ± 

u_ := POL¥_GCD(u, DIF(u, x>> 
If NUMBER? (u/u_ A 2) 

T(u/u_ A 2) u_ 


DO(u> := 

Prog 

If u = sqr(P) 

Prog 

PI := sqr(Pl) 

P2 := sqr(P2) 

o := ADJOIN<[± sfP], o> 
false 
Prog 

PI := SUBST (u, P, PI) 

P2 := SUBST (u, P, P2> 

o := ADJOIN< [u] , o) 

1 := ADJ0IN(P1, 1) 
r := ADJ0IN(P2, r) 

|"IMP IMP. IP-IP P2 ] 


UNDO(u) := 

Prog 

0 := REST (o) 

1 := REST(l) 
r := REST (r) 

PI := FIRST(l) 
P2 := FIRST (r) 


#5: play := [REUERSE(l), UECTOR{=„ k_„ DIM(l)), REUERSE(r), APPEND(REUERSE(o ) , [])]' 


A first example: 


J 2 x + 1 

#6: RECORD ( (3 - 7 = 20)) 

2 x + 1 

#7: 3 - 7 = 20 




i 


#13 : DO f— 1 = 

l 

#8: 

DO(P + 7) = 

2 x + 1 

3 = 27 


l 2 J 

X = 1 


#9: DO(LOG(P, 3)) = 


#10: DO 


#11: UNDO() = 


2 x +1=3 

1 

2 x +1 3 

2 2 


2 x +1=3 


#12: DO(P - 1) = 


2 x = 2 


#14: play = 


2 x + 1 

3 - 7 = 20 

2 x + 1 

3 =27 


2 x + 1 
2 x 


[P + 7] 
LN(P) 
LN{3) 
[P - 1] 
P 
2 
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Now let's see solving a quadratic equation step by step. (Don't forget the quote operator ‘ 
which is of essential importance. Josef) 


j -i 


#15: 

RECORD ( (x - 

- 4 x + 7 = 10)) 

#16: 

2 

x - 4-x + 7 

= 10 



4 

#17: 

DO(G - 3) = 

2 

x - 4 x + 4 

#18: 

D0(41?) = 

* 1 


|x - 2| = J7 


That was the "standard square root", let's try sqr 


#19: UNDO() = 

‘ x - 4 -x + 4 = 7 

#20: DO(sqr(P)) = 


1 1 

± sT? J 


#21: DO(P + 2) 


#22: play 


x - 2 = 

1 

x = ±47 + 2 

4 x + 7 


x - 4 x + 4 
x - 2 
x 


10 

7 

± 47 
±47+2 


[P - 3] 

[± m 

[0 + 2] 


I tried to add the features EXPAND and FACTOR. Doing so I had to change Johann ' s syntax: 

As my function is called doo, I use doo ( @ , "e" ) , doo ( 0 , "f " ) and doo ( 0 , "r " ) for expand- 
ing, factoring and drawing the square root. You are friendly invited to test your programming skills. 
I could not ’forbid” Derive simplifying the given equation to present unchanged in the first row of the 
play-table. Josef 


RECORD 

#46: 


b+x 4 b- (3 b-x) b-x V 


b-x 2 2 b+x 

b -2b x+x 

b+x 

4 

b (3 b-x) b-x 

#47: b-x 

b* 

! 2 b+x 

-2 b x+x 


#48: doo(G (b-x) (b+x))= 


2 2 3 

-(x+b) (x -4b x+11 b ) = (b-x) 


#49: doo(G,e)= 


3 22 3 3 223 

-x +3 -b-x -7b x-11 b = -x +3 -b-x -3b x+b 


#50: 

#51: 

#52: 

#53: 


3 2 2 

doo(G+x -3 -b-x +7-b -x) = 


-11 b 


doo(G-b )= 
G 

doo 


4 b 


3 2 

-12 b = 4b x 


-3 b 


iIEH 


2 3 

4 b x+b 


# 54 : 


2 2 
x -4b x+11 b 


2 

<x-b> 

2 2 
-(x+b) (x -4b x+11 b ) 

3 2? 3 

-x +3 -b-x -7b x-11 b 

3 

-11 b 

3 

-12 b 


b-x 


x+b 


3 

(b-x) 

3 2 2 3 

-x +3 b - x -3 b x+b 

2 3 

4 b x+b 

2 

4 -b -x 


[0 <x*b> (x-b) 2 ] 


[expand(G)] 

[ 3 2 2 

x -3 -b-x +7 b x+0 

R h 3 ] 

G 


2 

4 b 


-3 b 


x 
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Derive 6 acts as a teacher! 


Guiseppe Ornaqhi, Italy 

Using Derive 5.06 is it possible to simplify the following integral: 
INT(x A 3/(exp(x)-1 ),x,0, infinity) to (pi) A 4/1 5. 

But the function x A 3/(exp(x)-1 ) does not have closed form anti-derivative. 

I'd be interested to see how Derive works to find the exact form of the integral. 

Best regards 

Giuseppe 


Mike Law 

Maybe a couple of ideas (I am reaching ... these days I use tokens to ride the bus, I make many mis- 
takes counting change, tokens minimize my embarrassment) 

division by zero ... 

First guess is some sort of expansion ... 

I seem to recall Euler or Bernoulli numbers (exp(x)-1) A (-1) expansion of this 

if x/(exp(x)-1) mult, by x, then divide by x so that 1=(B0+B1x/1!+B2x A 2/2!+...)(1+x/2!+x A 2/3!+...) 
if x A 3 /() try mult. x A 3 then divide by x A 3 etc 

Not sure this is useful, however, ... a little more down memory lane I recall a topic in contour integra- 
tion 

Cauchy's theorem 

In the complex plane ... f(z) = 1/2*pi*i Integral on close path f(x)dx/(x-z) 

integral on closed path f(z) dz =2*pi*i SUM residues 

from: 

* Mathematical Methods of Physics by Mathews/Walker 1970 

They had some examples, none exactly like your problem, it is late at night and this got me out of bed 
when I felt something rising from the muddy bottom of the river of my past ... maybe I was dreaming. 

Regards, 

mike 


Albert Rich, 

Hello Giuseppe, 

The new Display Steps feature in the recently released Derive 6 makes it easy to answer your ques- 
tion. The definite integration rule applied to your integral is: 

If a>0 and p>-1, INT(x A p/(c*#e A (a*x)-c),x,0,inf) ==> p!*ZETA(p+1)/a A (p+1)/c 

Information about Derive 6 is available at www.derive.com and www.derive-europe.com . 

Aloha, 

Albert D. Rich 
Co-author of Derive 

DA/L: 


Follow Derive 6 stepwise to enhance your mathematical knowledge. You will find some of Mike Law's 
ideas realized. So Derive 6 acts as a teacher in mathematics. Johann Wiesenbauer provides some 
more information on this subject. 
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w 



Now I repeat the calculation - 
- starting stepwise simplification! 

o 

e * 1 

If a>0 and p >-I, 

p!-S(P + 1) 

CO 

P P + 1 

x a 

dx => 

a-x c 

J c-e - c 
0 

#3: 3!.£(4) 


If n>=0 is an integer, 
n ! 1*2*3* . . . *n 

# 4 : 6 -£( 4 ) 

If n>0 is even, 

n - 1 n 

2 • 7T 

£(n) => | BERNOULLI(n) | 

n! 

4 

4-6-tt 

# 5 : 

15 - 4 ! 

If n>=0 is an integer, 
n ! => 1*2*3* . . . *n 



Johann Wiesenbauer 


Well, as you can see from Albert's answer, Derive 
doesn't actually "compute" this integral, but 
essentially looks up the underlying formula in its 
memory. Just in case, you wonder where this 
formula comes from, here are some more lines, 
which might help to this end. (Note that all equa- 
tions below are the result of a simplification apart 
from the second one.) 

I hope these lines are fairly self-explanatory, but if 
this is not the case, please just tell me so and I 
will fill in the details. 

Cheers, 

Johann 


n :e Integer [0, ®) 
k :e Integer [0, ®) 
x : E Real (0, ®) 

n 

n ® - (k + 1) x x 

x - £ e = 

k=0 x 

e - 1 


CO 

‘ n ® 

x n®P A -(k+l)-x 

dx = x - £ J e dx 

x k=0 0 

e - 1 

0 


CO 

P n - (k + 1) -x -n-1 

Jxe dx=(k+l) -nt 

0 

to -n-1 

£ (k + 1) nt = nt £{n + 1) 

k=0 


Guiseppe 

I want to thank everyone who replied to my query. 

Mike Law 

It was a fun question, as usual I got more than I gave as your question took me back from a side trip 
to a wonderful book by Julian Havil called: "Gamma ... Exploring Euler's Constant" ... a very good 
book with plenty of detailed examples. 


Regards, 
mike law 




D 


-N-L#53 


Floating Point Error - A look into Derive ' s Interior 


P 23 


Floating Point Error? 


Valerui Anisiu, Romania anisiu@MATH. ubbcluj.ro 

Hello Derivers, 

I have found an annoying floating point error in Derive at least when working in high precision (it 
seems to be very old). 

E.g. if n:=10 A 50-1, then APPROX((n+10 A 23)/n,p) is simplified to 1 for p<=50, so that almost half of the 
digits are wrong! 

Using precision 51 one obtains a correct result, but changing APPROX((n+10 A 20)/n,51), 1 appears 
again! 

While writing this email I discovered that even for the implicit precision 10, 10000000500/9999999999 
and 20000001/20000000 are approximated to 1 instead of 1.000000050. 

It is very strange that the bug (feature?) was not discovered earlier! 

I found this error trying to compute with high precision an integral. Derive cannot approximate any 
integral with 30 digits in a reasonable amount of time (the Simpson's method is not suitable for high 
precision); even APPROX(INT(x A 4,x,0,1),30) is without hope [if you need it :-)]. So I have programmed 
a Gauss-like quadrature which was fast enough up to the precision 150 but less than half of the digits 
were correct due to the preceding FP error. 

Best regards, 

Valeriu 


Johann Wiesenbauer 

Hello Valeriu, 

As we all know, Derive has got a lot of strengths for which we love it, but it has also got some weak- 
nesses. I hate to say it, but the bad performance of APPROX() you reported in your posting (and I 
found myself many more examples) belongs definitely to the latter. 

I tried many different remedies, but none of them seems to work in all cases. In my opinion, for the 
time being the most promising one is to use a "self-made" approxQ like 

approx(e, p := 10) := INSERT < . , STRING(FL00R(e 10**)), -p) 

using the very powerful FLOOR(), which you can really rely on unlike APPROXQ. For example, you 
get 


50 

n := 10 - 1 


approx 


23 

n + 10 


n 


27 


1.000000000000000000000000001 


In the very crude form above, it works only for positive expressions e and displays trailing zeros, but I 
hope you got my point anyway. Maybe the underlying idea could also be used by Albert to resolve the 
current problems with APPROXQ in a more "professional" way. 


Cheers, 

Johann 
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Hello Derive users, 

Recently there have been several posting on the Derive Newsgroup email list concerning the accuracy of results 
returned by Derive. I felt the issue was of sufficient concern and interest that I wrote the following explanation: 

Rational vs. Floating Point Arithmetic 

Although Derive can display numbers in rational, decimal or scientific notation, internally all numbers are 
stored as integers or the ratio of two integers. There are numerous advantages, especially for a computer algebra 
system, to store numbers as rationals and use rational arithmetic, instead of storing them as floating point and 
using floating point arithmetic. These include: 

1 . Rational arithmetic is fundamental to mathematics and taught in elementary school. Floating point is a rela- 
tively recent invention designed for the convenience of calculators and computers. 

2. Floating point is unable to represent all rational numbers exactly (for example, 0.3333333 is not quite 1/3). 
The ability to represent ratios exactly is critical when symbolically simplifying expressions. 

3. The elementary operations (addition, subtract, multiplication, and division) can be performed exactly with- 
out roundoff errors using rational arithmetic. For example, the result of inverting a matrix of numbers twice 
using exact rational arithmetic is the original matrix. Whereas, successive floating point operations often result 
in the accumulation of roundoff errors. It is difficult to explain a result of 4.99999999 to students. 

Of course, there are disadvantages to rational arithmetic. These include: 

1 . Slower performance since arithmetic operations must be done in software rather than hardware (that is, until 
chip maker recognize the advantages of rational arithmetic and produce "floating slash processors"...). How- 
ever, note that floating point processors only provide arithmetic accurate to a fixed precision, but applications 
such as computer algebra require adjustable precision arithmetic. So the performance advantage of floating 
point is lost since the adjustable precision floating point arithmetic must be implemented in software. 

2. The more serious criticism of rational arithmetic is the phenomenon pointed out by Valeriu Anisiu and oth- 
ers. Unfortunately, the problem is not a bug in Derive but inherent in the uneven distribution of rational num- 
bers along the real number line. 

D. W. Matula and P. Kornerup present an algorithm for approximating rational numbers in their article "Finite 
Precision Rational Arithmetic: Slash Number Systems" published in the January 1985 issue of the IEEE Trans- 
actions of Computers. This algorithm, extended for adjustable precision arithmetic, was implemented in assem- 
bly language for muLISP and inherited by Derive ( muLISP is the LISP language interpreter in which Derive is 
written). 

Around small integers and simple fractions (for example, 1/2 and 2/3) there are relatively large intervals in 
which there is no better approximation for rational numbers than the integers or simple fractions themselves. 
Unfortunately rational numbers near the ends of these intervals may not be within the current digits of precision 
of the integer or simple fraction. Perhaps Johann Wiesenbauer can better explain the sparsity of rational num- 
bers whose numerators and denominators do not exceed a given size around integers and simple fractions. 

In the worse case around the integer 1 , rational numbers of only half the current digits of precision can be repre- 
sented. Therefore, the bottom line is: If you want to ensure that rational numbers of at least n digits are repre- 
sentable, use the Simplification tab of Derive ' s Options > Mode Settings command to set the digits of precision 
to 2*n. 

Of course, there is still the problem of catastrophic cancellation that occurs when approximating the difference 
of two nearly equal numbers (for details, see the description in the on-line help for the Precision Field of the 
Simplification tab of the Options > Mode Settings command). 

I hope this explanation and work-around helps. 

Aloha, Albert D. Rich, Applied Logician, Co-author of Derive 
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Hello Albert, 

First of all, many thanks for your taking the trouble to explain the internal workings of Derive, 
when it is trying to approximate an expression, which was certainly very instructive for many 
users. The bottomline of you message that you should use double presision in certain situa- 
tions is a little bit unsatisfactory though. 

Hoping not to arouse your anger even more, I have still some more questions. In the first 
place, you could use Maple inputting 

n : =10 A 50 ; 

evalf ( (n+10 A 23)/n,28) ; 
or Mathematica inputting 


n:=10 A 50 

N [ (n+10 A 23) /n,28] 

to get the correct result. I don't know, how they do this. Maybe they use a completely differ- 
ent approach or they use the same approach as Derive, but found a way to circumvent the 
problem. At any rate, this failure is not inherent to all CAS. (I know you didn't claim this, but 
the fact remains that Derive seems to be an exception here.) 

I remember that some years ago there was a quite similar problem with FLOOR() in Derive 
and I'm still very proud of the fact that pointing out this problem to you along with some sug- 
gestions how to fix it led to an extremely reliable implementation of FLOORQ. I have used it 
on many occasions since then without the slightest problems and - just to make up for the 
comparisons with Maple and Mathematica above - its performance is far superior to the per- 
formance of those rival programs on that score. 

I may be mistaken, but I think the situations with FLOOR and APPROX are quite similar in 
the way that certain numbers have to be dealt with more carefully. In the case, of FLOOR() 
these are the numbers very close to an integer. Please be patient with me and maybe this is 
completely unrealistic for some reason, but what about Derive detecting those cases on its 
own and taking counter measures like boosting the internal precision rather than leaving this 
to the user? 

There is one more thing concerning APPROX() that is "inconsistent" this time as I would call 
it. If you call APPROX from outside a program, you will always get a floating point number, 
which is very meaningful, indeed, and just what you expected. If APPROX() is returned by a 
program, it will be always a rational number represented as a fraction when using the default 
settings. I'm almost sure you will not understand me, but this tiny inconsistency is driving me 
crazy! There is simply no way of getting a floating point result other than adjusting the default 
settings which is committing sort of a sacrilege to me! 
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Albert's reply: 

Hello Johann, 

Thank you for your comments concerning this difficult problem. Requiring the user to use double 
precision to prevent this worst case error is certainly not optimal, but unfortunately I do not know how 
to fix the problem other than to internally use twice as much precision as the user specifies. I feel that 
would unnecessarily slow down the system. 

And yes, I did not claim that the problem is inherent to cas; however, I do claim that it is inherent to 
approximate rational arithmetic. I believe Maple and Mathematica use infinite precision rational arith- 
metic for exact calculuations, and adjustable precision floating point arithmetic for approximate ones. 
Derive uses infinite precision rational arithmetic for exact calculuations, and adjustable precision 
rational (aka floating slash) arithmetic for approximate ones. In addition to the considerable difficulty 
of implementing infinite precision floating point arithmetic in Pentium assembly language, there are 
innumerable problems (suffered by Derived competitors) caused by the need to convert from numbers 
from rational to floating point, and back. 

The APPROX problem is much more fundamental than the FLOOR problem. It is not just a matter of 
getting the right rational approximation, there is NO rational number that approximates the exact an- 
swer to within the current digits of precision whose numerator and denominator is small enough. 
Your idea of increasing the number of digits for the problem cases is intriguing, but knowing when 
and how much to increase the precision is beyond me. Perhaps if the result of approximation is an 
integer or simple fraction, and the approximation error exceeds the normal amount for the current pre- 
cision setting, then the precision could be increased. Sounds like a research project for a numerical 
analysis course. 

As I said in my previous email, internally all numbers in Derive are represented as rationals. There- 
fore, APPROX always returns a rational number. However, when an expression is simplified, a special 
check is made to see if APPROX is being called from the top-level. If so, the rational result it returns 
is displayed in decimal or scientific notation. For example, try simplifying [APPROX (pi) ] versus 
APPROX (pi) . 

As far as programming is concerned, the APPROX function is only able to set the precision to ap- 
proximate mode and the notation to scientific within the scope of its call. Therefore, if you want to 
approximate something and then analyze its decimal representation, include the analysis within the 
scope of the APPROX function. For example, APPROX (STRING (pi) ) simplifies to 
”3.141592653”; whereas, STRING (APPROX (pi) ) simplifies to ”1146408/364913”. 
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Josef, 


I noted with interest the problem of "Who makes himself a present?" on page 21 of issue #46. I re- 
cently encounter the same type of problem and also did a simulation to find the answer. But there is 
an exact solution to the matching problem as presented in Introduction to Probability Theory, by Hoel, 
Port, and Stone. 


As a reminder (from DNL#46): 

Who makes himself a present? 

At the occasion of a celebration 22 people bring a gift. All the gifts are distributed ran- 
domly among the 22 guests. Each of the participants receives one gift, so it can happen 
that one or the other will get back his own package. What is the average value of people 
getting back their own gift? 


The probability of no matches for n objects is: 

k 


ttl : P(n) 


n (-1) 
E 

k=0 kt 


And the probability of exactly r matches among n randomly permuted objects is: 

P(n - r) 

tt2 : Matches (r, n) := 

rt 

The probabilities of 0 to 22 matches among the 22 people bringing gifts is therefore: 

#3: t := TABLE(Matches(r, 22), r, 0, 22)' 


0 


#4: 


0.3678794 0.3678794 0.1839397 0.06131324 0.01532831 0.003065662 

6 7 8 9 10 

-5 -6 -6 -7 

0.0005109436 7.299195 10 9.123994 10 1.013777 10 1.013777 10 


11 


12 


13 


14 


15 


-9 -10 -11 -12 -13 

9.216155 10 7.680130 10 5.907788 10 4.219880 10 2.813063 10 


16 


17 


18 


19 


20 


-14 -15 -17 -18 -19 

1.759113 10 1.030867 10 5.857202 10 2.740211 10 2.055158 10 


21 


22 


-22 


0 8.896791 10 

Taking the dot product of the rows of the matrix gives the expected number of matches, which is 1 . 

t t =1 


# 8 : 12 

t is interesting to note that there are 0 matches for 21 . This is because if there are 21 matches, the 
last one is automatically a match also! 


But, what kind of probability distribution does this represent? If we take the limit of Matches() as n 
approaches infinity we get: 

-i 

e 

#9: lin Matches(r, n) = 

n -*to r! 

And this is simply the Poisson distribution with a mean of 1. 

Sincerely, 

Don Phillips 

phillipsm@gao.gov 

donphillips@starpower.net 
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Some comments on the change of publication form of the DERIVE Newsletter 


Lieber Josef, 

zum Ende der Print-Ausgabe der DNL mochte ich Dir und Deiner Gattin recht herzlich fur die vielen 
Stunden der Arbeit danken, die ihr beide in diese Publikationen gesteckt habt. Alle DNL waren stets 
liebevoll zusammengestellt und immer sehr interessant, so daft ich samtliche Exemplare aufbewahrt 
habe. 

Sicherlich werden auch die zukunftigen Ausgaben in gleicher Qualitat das Internet bereichern, worauf 
ich mich schon jetzt freue. 

Es erleichtert die Arbeit ungemein, wenn man auf die schnelle Verfugbarkeit von Artikeln zuruckgrei- 
fen und bei Fortbildungsveranstaltung darauf "verlinken" kann. Ich wurde es also sehr begruften, wenn 
auch altere Ausgaben so nach und nach in elektronischer Form erscheinen wurden, wenn ich auch 
die diesbezugliche Arbeit als nicht gering ansehe. 

Gregor Noll 


Dear Josef, 

I was reading the last issue of the DUG bulletin and I was feeling sorry that it will not be printed 
again. I know there are reasons that probably make this a wise decision, but receiving regularly the 
bulletin made me remember, from the distance, my very good friends in the first house of a little street 
in a small village in Austria... 

The lyrics of a Spanish song begins saying " Algo se muere en el alma, cuando un amigo se va... " 
("Something dies in your soul, when a friend goes... "). The DUG Bulletin was our friend and the con- 
nection between many friends for a long time. Now its somehow gone, but its spirit (and virtual .PDF 
version) will stay with us. 

Thank you for your (Noor's and Josefs) hard work along all these years. 

Your friend, 

Eugenio 


Lieber Josef, 

nachtraglich noch alles Gute zum neuen Jahr. Heute ist dein ’XXL’ Newsletter angekommen - immer 
ein hochsterfreulicher Posteingang, verspricht er doch stets interessante Lekture, viele Anregungen 
und neue Erkenntnisse. Ich kann allerdings nicht verhehlen, dass mich die Nachricht, den DNL zu- 
kunftig in elektronischer Form zu versenden, mich nicht besonders erfreut hat, wenngleich ich die 
Argumente naturlich nachvollziehen kann. Ich, fur meinen Teil, hatte lieber hohere Mitgliedsbeitrage in 
Kauf genommen, da ich nach wie vor gerne blattere. Deshalb habe ich die DNL ja auch, zu je drei 
Jahrgangen zusammengefasst, binden lassen. In diesen Banden kann man wunderbar herumblattern 
und entdeckt dabei immer wieder neue Oder in Vergessenheit geratene Dinge. Ich werde mir dann 
wohl den DNL immer selber ausdrucken mussen. 


Peter Luke-Rosendal 
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Dear Josef, 

I just received the latest DNL in the mail today. Thank you so much for the nice words 
regarding the Derive team in your "Letter of the Editor". We really enjoy working with 
you over the years! Derive wouldn't be the same product without all your questions, 
comments and suggestions. 

I too look forward to the new electronic era of the DUG, but backward with nostalgia at 
the printed versions. When you have new electronic information regarding the DUG 
membership and or a new DUG link, please let me know so I can update TI's Derive site 
at: http://education.ti.com/us/product/software/derive/features/groups.html . 

Additionally, many thanks for the article, "Wonderful World of DERIVE 6". You do a great 
job of explaining the new features and their application for classroom use. 

Aloha, 

Theresa 


Lieber Josef, 

danke fur das neue Heft. Die Ankundigung, es kunftig ohne Papier zu verschicken schmeckt mir sehr 
gut. Denn meine Ordner nehmen einfach nix mehr auf. Und die Hoffnung auf ein richtiges Register 
(und evtl. auch Links auf die *.dfw files etc.?) macht richtig geizig (d.h. geil, weil dieserja It. Werbung 
so ist). 

Allerdings muss ich auch eine kleine Riige mitschicken: 

Hast du bei der Danksagung im Editorial nicht die Noor vergessen? 

Deshalb besonderen Grull auch an sie. 

Viele Gruesse, 

Wolfgang 


Dear Josef, 

I am glad that the DUG membership is now free. After devaluation of our currency it has be- 
come hard to buy goods abroad. Please, take note of my E-mail so that you can send me the publica- 
tion dates of future DNL 's. 

Best regards, 

Marcelo. 


Dear Josef, 

Many thanks Josef for your work on the DNL, 

Please add me to the list for DNL publication date notification. 

I’ll miss looking forward to the orange envelope dropping through the door. 

On the other hand an electronic copy will lend itself to more rapid dissemination of information. 
I look forward to the next issue as always. 

All the best, 

Ian Pigram 
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Kaprekar's Sequence and his "Selfnumbers" 

Richard Schorn, Kaufbeuren (eMail richard.schorn@t-online.de) 

The "special number" 6174 (D-N-L #30, p 21 was discovered by D. R. Kaprekar (born in 1905), B.Sc. , 
S.T.C. of Devlali (India) in 1946. The introduction to his booklet "The New Constant 6174" reads as 
follows: 

61 74 is an invariant. It comes again and again from many four digital numbers by the 

"Reversed subtraction process " 

It is a matter of great pleasure and joy to bring this 61 74 from the working of the process. 

This is a new mathematical constant, (first discovered in 1946 by the author.) 

311 Devlali Camp 1-6-1959 D.R. Kaprekar 

Later on he writes: 

An article on this subject was first published by me in Scripta Mathematica of New York, vol XV, 
page 244-245. ...Iam the inventor of this new constant 6174 and ... 

Kaprekar did a lot of "number-crunching" (without computer or DERIVE!): 

...Hundreds of papers and note books were filled with several sorts of calculations and at last the 
secrets fully investigated. ... 

I learned about him by Martin Gardner's famous column "Mathematical Games" in "Scientific Ameri- 
can" in 1980 and I therefore bought several booklets from Mr. Kaprekar. 

The sequence 1, 2, 4, 8, 16, 23 (yes. 23, not 32!), 28, 38, 49, 62, ... seems to be Kaprekar's invention 
too: 

k 0 = 1 , k n = k n .i + sum of the digits of k n .i 

In his booklet "The Mathematics of the New Self Numbers" (1963) Kaprekar studied the problem 
whether a certain number appears in the sequence and how to find its predecessor. 3 has no genera- 
tor. He now (as an example) asks: 

Is there any generator for 86? The answer is No It will be seen by several trials that there is 

no generator for 86. 86 is the Self-Number. It is really Self-born. It is a <follows a Sanskrit ex- 
press ion> ...A number which has no generator is called a Self-Number. ... Thus between one and 
100 here are in all 13 Self-Numbers. They are 1, 3, 5, 7, 9, 20, 31,42,53, 64, 75, 86 and 97. ... 

There are contradictory opinions of the mathematical significance concerning Kaprekar's research, but 
in any case the results are well known in the world of "recreational mathematics" and give nice exam- 
ples for practicing computer programming. 

The TI-92-Programs use German words for the variables and inputs, but they are easily translated to 
other languages. 

Program KAPREKAR yields the sequence for different star- 
ting values. Program SELFNUMB gives the so-called Self- 
Numbers up to a predetermined value. Both programs need 
the function QUER (from german Quersumme) which calcula- 
tes the total of the digits of a given number. 

Here are some programs and their results to manage Kaprekars sequence. First one has to define a 
function which calculates the sum of the digits for a given number. Using a “normal” programming 
language one has to construct a loop using functions MOD and DIV (if available!). With DERIVE a 
very simple sum-of-the-digits-function is possible, following an idea of Josef Bohm in DNL #44, page 
33: 
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Q (n) := S (VECTOR ( CODES_TO_NAME (k) , k, NAME_TO_CODES ( n ) ) ) 

A step-by-step demonstration with the number 471 10815 will clarify the process: 

• NAME_TO_CODES (47110815) = [52, 55, 49, 49, 48, 56, 49, 53] 

One gets a vector with the ASCII-Codes of the eight digits. 

• VECTOR (CODES_TO_NAME (k) ,k, [52,55,49,49,48,56,49,53]) = [4 , 7 , 1 , 1 , 0 , 8 , 1 , 5] 

The “ASCII-names” of the eight digits make up a new vector. 

• SUM ( [4 , 7, 1, 1, 0, 8, 1, 5]) =27 

The sum (Z) of the elements is determined. 

N (n ) := n + S (VECTOR (CODES_TO_NAME (k) , k, NAME_TO_CODES (n) ) ) 

Obviously this function determines the next element in Kaprekars sequence. Without the help of com- 
puters Kaprekar discovered remarkable numbers with three and four generators: 

N( 999999999999999999999893) = 1000000000000000000000102 
N( 999999999999999999999902) = 1000000000000000000000102 
N(1000000000000000000000091) = 1000000000000000000000102 
N(1000000000000000000000100) = 1000000000000000000000102 

In [3] he writes on page 18: 

This number will have a historical value in future as discovered by D. R. Kaprekar 
on 7-6-61 by his methods and also discovered by Prof. Gunjikar on 7-9-61 by his 
different methods. 

With function n (n) using iterates one gets the sequence up to a given index: 

ITERATES (N (n) , n, 1, 25) = [l, 2, 4, 8, 16, 23, 28, 38, 49, 62, 70, 77, 91, 
101, 103, 107, 115, 122, 127, 137, 148, 161, 169, 185, 199, 218] 

In a similar way we calculate the m th element of the sequence: 

KAPREKAR (m) := ITERATE (N (n) , n, 1, m) 

Some examples: 

KAPREKAR (10) =70 
KAPREKAR (100) = 1213 
KAPREKAR (10000) = 213703 
KAPREKAR (99999) = 2609882 
KAPREKAR (100000) = 2609917 

The function self ( n ) returns true, if the argument n has no predecessor or in Kaprekars words if n 
is a “self-number”. 

SELF (20) = true 
SELF ( 70 ) = false 
SELF ( 1000000 ) = true 
SELF (2609917 ) = false 


The last program yields a list of selfnumbers up to a given number: 


SELFNUMBER ( 

1000) 

= [1, 

3, 

5, 7, 

9, 2 C 

), 31, 

42, 

53, 

64, 11 

5 , 86 , 

97, 

108 , 

110 , 

121 , 

132 , 

143, 

154 , 

165, 

176, 

187 , 

198 , 

209, 

211 , 

222 , 

233, 

244 , 

255, 

266, 

277 , 

288 , 

299, 

310 , 

312 , 

323, 

334 , 

345, 

356, 

367 , 

378 , 

389, 

400 , 

411, 

413, 

424 , 

435, 

446, 

457 , 

468 , 

479, 

490 , 

501 , 

512 , 

514 , 

525, 

536, 

547 , 

558 , 

569, 

580 , 

591 , 

602 , 

613, 

615, 

626, 

637 , 

648 , 

659, 

670 , 

681 , 

692 , 

703, 

714 , 

716, 

727 , 

738 , 

749, 

760 , 

771 , 

782 , 

793, 

804 , 

815, 

817 , 

828 , 

839, 

850 , 

861 , 

872 , 

883, 

894 , 

905, 

916, 

918 , 

929, 

940 , 

951 , 

962 , 

973, 

984 , 

995] 












Note the the curious part of this list, starting with 21 1 ! 
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TI-92 -Programs KAPREKAR. 9XP and SELFNUMB . 9XP using Function QUER.9XF 


quer (x) 

Func 

Local h,s 

0+s 

x^h 

While h^0 

s+remai n (h , 10)-»s 
intDiv(h,10)->h 
EndWhile 
s 

EndFunc 


(t 1 1 gebraTca fclotherlpr gw I ole 1 ear. Up 


D 


■ guer( 1 23005689732589) 


|guer<123QQ5689732589> 


68 


kaprekar ( ) 

Prgm 

Local aus,z,g 
ClrlO 
0->z 
" 'Aaus 

Input "Initial Number", g 
ClrlO 

While z<79 

aus&" "&stri ng (g)^aus 
If dim(aus)>34 Then 
Disp aus 
" "-»aus 
Endlf 

g+quer (g)->g 
z+1+z 
EndWhi 1 e 
Disp aus 
EndPrgm 





IA;V 

3? 

“T FE 

■;-|PrgrAl0 

37 


n 

5 10 11 

13 

17 25 

32 

37 

47 58 71 

79 


95 

109 

19 

130 134 

142 

149 163 

173 


184 

197 

214 

221 

226 

236 

247 

260 

268 


284 

298 

317 

328 

341 

349 

365 

379 

398 


418 

431 

439 

455 

469 

488 

508 

521 

529 


545 

559 

578 

598 

620 

628 

644 

658 

677 


697 

719 

736 

752 

766 

785 

805 

818 

835 


MAIN 


RAD AUTD 



FUNC H/3D 




selfnumber() 

Prgm 

1 ocal q , f 1 ag , i , n 
{1}+selfn 
{1 }->fol ge 
l^sel fz 
ClrlO 

input "Upper Bound?", n 
For i , 2 , n 
1 -*-f 1 ag 

for k, 1 , sel fz 
if i>folge[k] 
fol - 

ge[k]+quer(folge[k] )^folge[k] 
if i=folge[k] 

0->f 1 ag 
EndFor 

If flag=1 then 
sel fz+1^sel fz 
augment (sel fn , {i } )^sel fn 
augment (fol ge , {i } )->fol ge 
endlf 
endFor 

ClrlO 

"D.R.Kaprekar's "&char(34)+tt 
tt&"Sel fnumbers"&char (34)^tt 
Disp tt 

disp "Upper Bound "&string(n) 
disp string(selfz)&" Sel fnumbers : " 
" "->aus 

for i , 1 , sel fz 

aus&" "Sstri ng (sel fn [i ] )->aus 
if mod(i ,8)=0 then 
Disp aus 
" "->aus 
Endlf 
EndFor 
Disp aus 

EndPrgm 


[1] D.R. Kaprekar, The New Constant 6174, Devlali, 1959 

[2] D.R. Kaprekar, The Mathematics of the New Self Numbers", 1963 
(This is the first booklet of a series with the same title) 

[3] D.R. Kaprekar, Self numbers and numbers with 3 and 4 generators, Devlali, 1962 
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Hallo Josef, 

als Dank fur deine schone Nr. 51 mochte ich dir eine kleine Knobelei schicken: 
As a thank you for fine DNL#51, I'd like to send a little brain twister. 


SEP - Zahlen SEP - Numbers 

Riideger Baumann, Celle 

Nimm eine natiirliche Zahl >10, zerlege sie in ihre Ziffem und bilde deren Summe. 
Take any natural number >10 and calculate its sum of digits 


Beispiel: 28 — > 2,8: Summe: 10. 

Nun hange die Summe hinten an und bilde wieder eine Summe aus zwei Summanden, indem du die 
erste Ziffer weglasst. 

Attach the sum and calculate again the sum of all numbers after removing the first one. 

Im Beispiel: 2, 8, 10, 

Auf diese Weise machst du weiter: 2, 8, 10, 18, 28 - oh Wunder: die Startzahl ist wieder erschienen. 
Dies ist natiirlich ein seltener Gliicksfall. 

Dir zu Ehren soli die 28 und jede solche Zahl SEP-Zahl (sich selbst produzierende Zahl) heiGen. 

So you go on: 2, 8, 10, 18, 28 - oh what a miracle, the initial number appears again. This is a rare 
lucky chance. To honor you, III 28 and all such numbers a SEP-number (self producing number). 

(That needs some explication: Sepp is an Austrian and Bavarian form for Josef, Guiseppe is its Italian 
origin.. In Vienna - and Eastern Austria we have another form for Josef: Pepi or Peperl.) 


Zweites Beispiel (langsam): Startzahl 1104 

1+1+0+4 = 6 
1 +0+4+6 =11 
0+4+6+ 1 1=21 
4 + 6+11+ 21 =42 
6+11+21+42 = 80 


11 +21 +42 + 80= 154 
21 +42 + 80+ 154 = 297 
42 + 80+ 154 + 297 = 573 
80 + 154 + 297 + 573 = ???? 1104 

BINGO!! SEP-Zahl /Number 


Herausforderung: Man schreibe ein DERIVE-Programm zur Ermittlung von SEP-Zahlen. Wie viele 
davon < 10 10 gibt es? 

Challenge: Write a DERIVE program to find SEP-numbers. How many of them can be found < 10 10 ? 


I tried to do my “home work”: 

sep(11436171) 

[1, 1, 4, 3, 6, 1, 7, 1, 24, 47, 93, 182, 361, 716, 1431, 2855, 5709, 11394, 22741, 
45389, 90596, 180831, 360946, 720461, 1438067, 2870425, 5729456, 11436171] 

5 6 

seppp(10 + 1, 10 ) 

10 

[120284, 129106, 147640, 156146, 174680, 183186, 298320, 355419, 694280, 925993] 
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Sierpinski-Tetraeder: Darstellung in 3D unter Derive 5 
AG CoMa, Gutenbergschule Wiesbaden 
2.Hj. 2001/02 

www.gutenberg-gym.de/coma/ 

1. Berechnung der Eckpunktskoordinaten des Tetraeders; k ist hierbei die Kantenlange 
Calculating the coordinates of the vertices of a tetrahedron with edge = k 


a := [0, 0, 0], b := [k, 0, 0], c := 

1 k 1 

k, sT3, 0 

ii 

■ ■ 


2 2 




Die Herleitung, durch Anwendung des Satzes von Pythagoras, zeigen wir mit der Kantenlange k = 5. 
Die Grundflache des Tetraeders ABC sieht so aus: 


Firstly we apply the Pythagorean Theorem and work with k = 5. We start with the base triangle ABC. 


#2: dreieck 


0 0 

5 0 

2.5 2 .5 - \T3 

0 0 


Die Spitze des Tetraeders D liegt uber dem Schnittpunkt der 
Seitenhalbierenden: 

The top of the tetrahedron D is above the intersection point of 
the perpendicular bisectors of the sides. 


seitenhalb 


[5, 0] + [2.5, 2 .5 - \T3] 

[ 0 , 0 ] 

2 

[2.5, 2 .5 - \T3] 

[5, 0] 

2 

[ 5 , 0 ] 

[2.5, 2.5 - \T3] 

2 

D1 (k/2, kV3/6) 





Das Dreieck AD^ ist rechtwinklig und es gilt: k 2 = dz 2 + 1 AD 1 1 2 = dz 2 + 

Wir losen nach dz (der z-Koordinate von D) auf. 

Triangle AD is a right triangle and we solve the above expression vor dz = 3rd coordinate of D. 


tt4: 
# 5 : 


#6 : 


In put Mode := Uord 


2 2 
k = dz 

f 2 

SOLUE k 



\T6 k \T6 k 

#7 : dz = - v dz = 

3 3 


Damit sind die Koordinaten von D bekannt: 


D 


fk kV3 kVeH 


So we know the coordinates of vertex D: 


J 
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k := 3 

Nun haben wir unser Tetraeder schon 3D darge- 
stellt, aber wir wollten es am Ursprung zentrieren. 

Now we see our tetrahedron in a 3D-view, but we 
would like to have it centered in the origin. 


4. Wir verschieben nun unsere Punkte a,b,c,d, so dass das Tetraeder zentriert liegt. 

We translate points a, b, c, d such that the origin will become the center. 

Dazu subtrahieren wir den Vektor von M von jedem Punkt. 

To do so we perform 4 substractions: point minus vector M 

tetraederl := UECTOR(UECTOR(teti*aeder — n, j, 1, 4), i, 1, 4) 
#22 : i,j 



Origin O is center of the tetrahedron 
with edge k = 3. 


Wir automatisieren jetzt die Darstellung des Tetraeders. 

Zur Vereinfachung nehmen wir hier fur die Kantenlange k=3 an. 

Die Definition des Tetraeders: 

We automate the presentation of the tetrahedron, to make it easier we take k = 3 and repeat the 
definition: 


#23: 1 

t := 3 

r i 

k 1 


#24: 

a := [0, 0, 0], b := [k, 0, 0], c := 

k, 

2 

vT3, 0 

2 

, d : = 


#25: [f 1 := [a, b, d, a], f 2 := [a, c, d, a], f 3 := [b, c, d, b], f 4 := [a, b, c, a]] 

#26: tetraeder := [fl, f 2 , f 3 , f4] 

Da wir haufiger verschieben wollen definieren wir uns eine Funktion trans () , die fur die Verschie- 
bung eines Objektes mit Vektor v, wobei das Objekt in der Flachendarstellung vorliegt, zustandig ist. 

We intend to apply a translation more often, so we develop a function trans (), which is responsible 
for translating an object given by its faces by a vector v. 
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trans(obj, u) := UECTOR(UECTOR(obj + u, j, 1, 4), i, 1, DIM(obj)) 
#27: i,j 

#28: trans(tetraeder, b - a) 

Jet zt wird das Tetraeder in drei Richtungen verschoben: 

Now we translate the tetrahedron into three directions: 


#30: [tetraeder, trans (tetraeder, b - a), trans (tetraeder, c - a), trans (tetraeder, d - 

a)] 


0 < x,y,z < 6 

A click with the right mouse key into each single 
tetrahedron opens the Edi t window and one 
can change the Plot Color using Scheme 
Custom. 



Fur die Verschiebung des aufzubauenden Sierpinski-Tetraeders braucht man eine modifizierte Ver- 
schiebungsfunktion, weil die Struktur des Sierpinski-Tetraeders eine Liste von Listen ist. 

For translating the requested Sierpinski-tetrahedron we need a modified translation function, because 
the structure of the Sierpinski-tetrahedron is a list of lists. 

transl(obj, u) := UECTOR(UECTOR(UECTOR(obj + u, k, 1, 4), j, 1, 4), i, 1, 

#31: i,j,k 

DIM(obj)) 

Nun endlich kann das Sierpinski-Tetraeder rekursiv aufgebaut werden. 

Zu Beginn (n=0) steht das Orginal-Tetraeder und seine 3 verschobenen kleinen Bruder. 

Dieses Objekt wird nun als Ganzes urn die doppelte Kantenlange in 3 Richtungen verschoben. 

Finally we can build the Sierpinski tetrahedron recursively. 

We start with n = 0 - the initial tetrahedron - and its three translated little brothers (or sisters). 

Then we translate this object as a whole by double side length into three directions 


sier(n) := 

If n = 0 

#31: [tetraeder, trans (tetraeder, b - a), 

APPEND{ [sier(n - 1), transl(sier(n - 


transftetraeder, c - a), trans(tetraeder, d - a)] 

1), 2^n(b - a)), transl (s ier(n - 1), 2 A n - (c - a)), transl (s ier(n - 1), 2 A n(d - 


a>)]> 
to edit as: 

sier(n) := IF(n = 0, 

[tetraeder, trans(tetraeder, b - a), trans(tetraeder, c - a), trans(tetraeder, d - a)], 
APPEND([sier(n - 1), transl (sier(n - 1), 2 A n (b - a)), transl (sier(n - 1), 2 A n (c - a)), 
transl (sier(n - 1), 2 A n (d - a))]) 


) 
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Teste die folgenden Beispiele / Test the following examples : 



Nun sollte das Tetraeder am besten noch den Mittelpunkt im Ursprung haben! 
Alles kein Problem, verschieben wir es also. Final translation into the origin. 


#37: sier_c(n) 


transl sier(n), - 2 



sT3 

* 

6 


T6 

12 


1) 


sier_c(3) 


sier_c(4) 




view from above 



view from the bottom 
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Titbits from Algebra and Number Theory (27) 

by Johann Wiesenbauer, Vienna 

In the current issue of the "Titbits" I'll try to keep an old promise I gave Peter 
Schofield back in December 2002. At that time he sent me a short email in which he 
was "wondering if DERIVE would be strengthened by including some kind of decimal no- 
tation corresponding to rational numbers". It soon became clear by the programs he 
included that he meant some kind of conversion of a rational number to a form which 
makes visible the preperiod and the period of it. I wrote back that I found his idea very 
interesting, but was still thinking about an optimal implementation. I particular, I was 
aiming at some sort of a "Swiss Army Knife", which could do all sorts of conversions in 
terms of changing to different bases. In fact, I was trying hard for some time to 
achieve this goal, but no solution looked really satisfactory to me then. Sadly enough, I 
forgot about this matter eventually. 

Now, as some of you might know, quite recently Volker Loose has brought up exactly the 
same issue on the Derive-newsgroup again. And once more I looked into this problem, 
but this time I succeeded in writing a routine that comes close enough to my picture of 
a "Swiss Army Knife" to be worth being published on the newsgroup. In the meantime, I 
have even added some more improvements (hopefully you agree that these are improve- 
ments, indeed!) and you can see the outcome converter, k,b) on the next page along with 
some examples. 

Basically, you can use this utility function to convert a positive rational number to an- 
other representation, where you can see both the preperiod and period of r w.r.t. any 
base b with the usual restriction 2<=b<=36. As a spin-off, you may also use this function 
to convert positive integers to other bases b than 10, also with a fixed block-length k, 
as this is often needed in coding theory. The parameters r, b and k have the following 
meaning: 

r is the nonnegative rational number to be converted, 

b is the base for all calculations and the output ( 2<=b<=36; b-10 by default), 
k is sort of a "switch". To be more precise, 

• if k and r are both positive integers, then the output consists of the k least sig- 
nificant digits of r 

• if k=0, then the output consists of as many digits as necessary, which is the de- 
fault case 

• if k=-l, then the output is only the period length of r 

• if k=-2, then the output is only the preperiod length of r 

• if k=-3, then the output is a vector containing both the preperiod length and the 
period length of r 
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converter, b ;= 10, k := 0, a_, d_, n_, r_ := 0, s_) := 
Prog 

If r = 0 

RETURN IF(k < 0, 0, REST(STRING(10 A k))) 
n_ := NUMERATOR(r) 
cL := DENOMINATOR(r) 

Loop 

a_ == GCD(d_, b) 

If a_ = 1 exit 
r_ : + 1 
n_ : * b/a_ 
d_ :/ a_ 

S_ := b 

Loop 

If M0D(s_, d_) < 1 exit 
s_ : * b 

n_ :* (s_ - l)/d_ 

OutputBase ;= b 

a_ == STRING (FLOOR (n_, s_ - 1)) 

Loop 

If DIM(a_) > r_ exit 
a_ := ADI0IN("0", a_) 
a_ == INSERTC" . " , a_, -r_) 
s_ == REST (STRING(s_ + M0D(n_, s_ - 1))) 
OutputBase ;= Decimal 
If k = -1 

RETURN DIM(s_) - 0 A (d_ - 1) 

If k = -2 
RETURN r_ 

If k = -3 

RETURN [r_, DIM(s_) - 0 A (d_ - 1)] 

If d_ = 1 
Prog 

If r_ > 0 
RETURN a_ 

a_ == DELETE (a_, -1) 

If k = 0 
RETURN a_ 

Loop 

If DIM(a_) = k 
RETURN a_ 

If DIM(a_) < k 

a_ == ADIOINC'O", a_) 
a_ == REST (a_) 

n_ := APPEND(a_, "<" , s_, ">") 


convert 


1 "1 


28 J 


0.03<571428> 


convert(13, 2) = 1101 


convert 


1 

, 2 

3 J 


0 .< 01 > 


p= 2231588810593399 is the smallest prime such 1/p has period 2003 (C f. 
www.primepuzzles.net and puzzle 208 there.) Indeed, we get: 
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#5: 

# 6 : 


convert 


2231588810593399 


0. <000000000000000448111226966625292420798398833515682905127413360- 
6093400509692366772061856346118381238015662344060034958096567671- 
8465239302279637019360464382481161390973304491412032169179317368- 
4363304603984145005324888834688688141604620765892780956425386302- 
7921045155758188142921476163065425554949408238482149451490093731- 
7548342921070719501851040433529102728286338723009466465331893113- 
8297471379355688192826936599518389164526410843653223163360019280- 
6511374821246283580552313564246559241032247330963586193127370501- 
8194689153289081804961825447591063042310608638405893039308877640- 
5753492360271462100989507408730431522985979959583881619285805717- 
4956985794545949343822990493183545432966827488090387733838229816- 
6874784898242373565302056595598059535989834556216101551511150280- 
1606335958025171245075306417378727690982587588301548934406714991- 
4916403128781830242186951419655355421202318965031001893431061122- 
5192071835842486485977183386982311206045204856823054908205090005- 
1457095252402440185349219597697282968516959889113713793339827398- 
0169042709280151940426635818586478183624967375739407827647070489- 
8731124298744735851336865455313228325285689895982135213079015639- 
3973517005569009365310608190240051620674573829012580756056830914- 
0861091665055521167919693410936946690301528609529890600515052971- 
4051434420936889742716746533408914756039524769207326344980705363- 
1062052234605625785530469775250398735345956785311378483544651453- 
6175180821393428602003951104683777930590963511671871505368152083- 
7999021520569561660741408127099234360647041167779653106470237958- 
5214299697314497545770214930786241140134761337903093696212628931- 
2954218500487228960441587275257453254772998479084037652756473151- 
9208332955045286516205950157354439047655931391572096101419557550- 
7988570568757731946625834990006145506154130277608644778693485112- 
5051414053848753593570942608521413148260947750152547980381398390- 
7679025224708477593900111958464815686165291730713241052811379472- 
9517308487175627422606797290336183884481746804928823750331514216- 
26685329883038153401> 


convert| , 10, -1 1 = 2003 

^ 2231588810593399 ) 


convert , 10, -21 = 0 

l 2231588810593399 J 

convert] , 10, - 3 ] = [0, 2003] 

l. 2231588810593399 ) 


As you might have noticed I haven't mention any computation times so far. The reason 
for this is very simple: It takes usually only fractions of a second to get those results. 
For example, the last 4 computations took 0.08s each on my 2GHZ-PC. Apart from all 
those "bells and whistles" in my program, which you might find useful or not, this may 
well be the biggest difference to Peter's sample programs mentioned above (dealing 
with the basic case b=10 only), which are programmed in a completely different way and 
are say “ a little bit slower", though giving the correct results. 

Let's close this topic with a look at the period lengths of 1/p for all primes p below 100. 
These are given by the following table: 
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f 1 


A 





A 


TABLE 

convert 

— 

, 10, 

-1 

. P. 

SELECT (PRIME(q_) , q_, 1, : 

100) 





1 P 


J 





) 


' 2 : 

3 5 

7 

11 

13 

17 

19 

23 

29 

31 37 41 43 47 

53 

59 

. o : 

L 0 

6 

2 

6 

16 

18 

22 

28 

15 3 5 21 46 

13 

58 

61 

67 

71 

73 

79 

83 

89 

97 ' 
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You will immediately see the interesting property, that the period length of 1/p is al- 
ways a divisor of p-1, if p doesn't divide 10. Only for a few primes p, whose values you 
can conclude from the table, the maximum p-1 is actually reached, though. What do 
those primes have in common? Well, it's not my intention to put you on tenterhooks: 
These are exactly the primes for which 10 is a primitive root, i.e. a generator in the 
cyclic prime residue class group mod p. 

SELECT (PRIMITIVE_R00T (p , 10) = 10, p, SELECT (PRIME(q_) , q_, 1,100)) 

[7, 17, 19, 23, 29, 47, 59, 61, 97] 

Hope to see you all in Montreal. And as always, if you have any comments or suggestions, 
please let me know! (i.wiesenbauer@tuwien.ac.at) 
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Hi fellow Derivers, 

I am hoping to set up an international group of mathematicians, who are now, or will be, interested in 
promoting the many mathematical achievements of Ernesto Cesaro. 

He seems to be known mainly for the Cesaro summation formula for certain divergent series. This 
single achievement, is only a minuscule 'drop in the bucket' of his vast host of accomplishments. 

The centenary of his death will be on 12 September 2006. As a group we could publish papers, that 
bring his work before the mathematical diaspora, not merely for the purpose of reviving his name, but 
also because of its innate worth. There is also a great deal of untapped potential in his writings, and 
we all can be beneficiaries. I am sure that many college curricula, at least, would benefit with some of 
his findings. 

Personally, I have found so much advantage from Intrinsic Geometry, that it seems to have no end. 
What I have in mind is that there could be occasional talks/lectures on campus to 'sound out 1 the lo- 
cals and we could enlist all interested parties in a Cesaro subgroup on a list server, either Jiscmail or 
Yahoo Groups. Eventually, two years down the track, we can have a centenary celebration, the nature 
of which is yet to be determined. 

I have included a short biography below. Please let me have feedback as to your own ideas and also 
please tell anyone that you think may be interested, even if they are not fellow Derivers. I welcome all 
suggestions. 

Many thanks, David Halprin (convenor of eDUG) 
http://groups.yahoo.com/group/eDug 


CESARO, Ernesto 

Mathematician, born in Naples on 12th. March 1859. At the age of 14, due to unforseen circum- 
stances, he was obliged to go and live with a brother at Lieges, where he attended the "School of 
Mines", and, because of his unusual approach to mathematical research, very soon attracted the at- 
tention of Catalan, who then became his close friend and protector. Not very drawn to studying engi- 
neering, he succeeded after many trials to begin the study of mathematics at the University of Rome, 
but because of poor health and financial difficulties was not able to attend lectures regularly, nor to 
graduate. However, even in this difficult period, he succeeded in publishing about 100 works, which 
helped him to gain, in 1886, the Chair of Analytical Algebra at the University of Palermo. In 1891 he 
went to Naples to teach Infinitesimal Calculus, and there, while on the point of taking up the Chair of 
Rational Mechanics at the University of Bologna, in an attempt to save his son from drowning, died 
with him, (12th. September 1906). 

An original and most prolific mind, which soared into most widely dispersed fields:- 

From elementary geometry via Intrinsic Coordinates, [including two conditions for immobility in a pla- 
ne, (of a point and of a line)] to the theory of numbers, etc. 

From analysis, both algebraic and infinitesimal, to the calculation of probability, u.s.w. 

From differential geometry to mathematical-physics, etc.. 

I have a set of three volumes of most of his papers "Opere Scelte", (Edizioni Cremonese, Roma 
1968), but unfortunately many of his best papers on Intrinsic Geometry are not therein, but are dis- 
persed throughout European Universities. (I have most of these too.) 

The partly self-taught nature of his scientific development, far from diverting his research activity, led 
him into every field; to ask and answer new questions, and cast new light upon problems, already very 
old. He was a definitive polymath (polyMATH), in every sense of the word. 

The work of Cesaro, because of its great volume, and the fragmentary nature found in every part, 
makes it difficult to draw a conclusion. It is enough to remember his contributions to asymptotic arith- 
metic, particularly "Excursions in Arithmetic of the Infinite" in Annals of Mathematics (Milan 1885), the 
"Course of Algebraic Analysis" (Turin 1884), which reflects the arithmetical view of algebra, character- 
istic of Cesaro, and which was then translated into German; the "Introduction to the Mathematical 
Theory of Elasticity (Turin 1894); the "Lessons in Intrinsic Geometry" (Naples 1896), a distinctly per- 
sonal work of Cesaro, and also this was translated into German, and finally the "Elements of Infini- 
tesimal Calculus" (Naples 1899, 2nd. edition in 1905). 






